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We report on a study of electron-hydrogen scattering, using a combination of a 
modified method of polarized orbitals and the optical potential formalism. The 
calculation is restricted to P waves in the elastic region, where the correlation functions 
are of Hylleraas type. It is found that the phase shifts are not significantly affected by the 
modification of the target function by a method similar to the method of polarized 
orbitals and they are close to the phase shifts calculated earlier by Bhatia [Phys. Rev. A 
69, 032714 (2004)]. This indicates that the correlation function is general enough to 
include the target distortion (polarization) in the presence of the incident electron. The 
important fact is that in the present calculation, to obtain similar results only 3 5 -term 
correlation function is needed in the wave function compared to the 220- term wave 
function required in the above-mentioned previous calculation. Results for the phase 
shifts, obtained in the present hybrid formalism, are rigorous lower bounds to the exact 
phase shifts. 

PACS numbers: 34.80.Bm 


I. Introduction 

Scattering of electrons by hydrogenic systems is always of interest because the 
target function is known exactly, allowing us to test the various methods of calculation. 
At low incident energies, the distortion of the target produced by the incident electron is 
important. One of the methods used to take into account this distortion is the method of 
polarized orbitals [1], which includes the effect of polarization and essential physics in 
the ansatz for the scattering wave function. Various other approximations have been 
used: Kohn-Feshbach variational method [2], Kohn variational method [3], R-matrix 
method [4], and the finite element method [5]. hi previous papers [6, 7], the P - wave e-H 
and e-He + scattering phase shifts were calculated by using the Feshbach projection 
operator formalism [8]. The results obtained are accurate and have variational lower 
bounds. But it has not been possible to take into account in the Feshbach formalism the 
distortion produced by the incident electron which results in a direct -1/r 4 potential in the 
scattering equation. 

In [9], the S-wavc electron-hydrogen scattering phase shifts were calculated by 
using a hybrid method in which both long-range potential proportional to -1/r 4 and short- 
range correlations via an optical potential were included in the scattering equation at the 
same time. We did not use any projection operators in this calculation [9] but the 
important property, namely, that the phase shifts are rigorous lower bounds to the exact 
phase shifts, is retained [10]. We follow the same procedure in the present calculation on 


P - wave scattering as in [9]. We use Rydberg units: energy in Rydbergs and length in 
Bohr radius a 0 . 


II. Theory 

The total spatial wave function for e-H partial wave (denoted by L) problem is 
written as 





N is the number of terms in the expansion. The target wave function is given by 


( 1 ) 


^ 0 (r 2 ) = (Z 3 /^)°- 5 e^ 


( 2 ) 


where Z is the charge of the nucleus. The ( ± ) above refers to singlet (upper sign) or 
triplet (lower sign) scattering, respectively. Beyond the terms containing u , explicitly 
(those are the terms giving rise to the exchange approximation), the functions ® L are the 
correlation functions. These functions include all the many-body effects and the 
scattering equation is a single-particle equation. For arbitrary L this function is most 
efficiently written in terms of the Euler angles [11]: 

&L = Ul + ( r \’ r 2’ r n) D L K+ + j l~ (ri,r 2 ,r u )D L K - {&,(/>,¥)\ ( 3 ) 

The D k e functions (s=+l,-l) are called rotational harmonics [11] and are functions of 
the Euler angles . The/s above are generalized “radial” functions which depend 

on the three residual coordinates that are required (beyond the Euler angles) to define the 
two vectors iq and r^. The distance between two electrons is r l2 =\r ] -r 2 | . The radial 
functions are defined as follows 


/i' + = cos (*9i2 / 2 )[/ (q , r 2 , r n ) ± f(r 2 , r x , r l2 )] , 
/' = sin(t9 12 / 2 )[/ (q ,r 2 ,r l2 ) + f(r 2 , r x , r l2 )] , 


9 X2 is the angle between r x and/y , and 


f(r l ,r 2 ,r X2 ) = Y J C t 
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The minimum value of l is equal to 1 while that of m and n is 0. First, we derive 
the wave function of the scattered electron without the long-range correlations and it is 
given implicitly by 


(V) 


\[Y; 0 {Q. x )f Q {r 2 ){H-E)V L W2=Q, 

H is the Hamiltonian and E is the total energy of e-H scattering system: 


2 2 2Z 2Z 2 

(8) 

h = -v 2 x -v 2 2 +— , 

r, r 2 r X2 

E = k 2 -Z 2 , 

(9) 


2 

k is the kinetic energy of the incident electron and Z= 1 is the nuclear charge which is 
equal tol in the present calculation (electron-hydrogen scattering). 

In order to derive the equation eq. (7) for the scattering function u(r\) = u L (r, ) , the 
coefficients C,_ must be known. Taking N= 1, for illustration, we calculate the functional 


(arising in the Kohn variational principle) 

/ =< W L (k x ,P 2 ) \H-E\ T z (F 15 F 2 ) >, (10) 

which can be written as 

I = A + C X B + C X D, (11) 

where 

A =< (f 2 )u(r x ) ± (1 ^ 2)] | H - E | [0 O (r 2 )u(r x ) + (1 2)] > , (1 2) 

and 

B = 2<<S>f{r l ,r 2 )\H-E\[<p Q (r 2 )u(r x ) ± (1 <-> 2)] >= 4 < V x (r, )u(r x ) > , (13) 

where 

VM) =< I H -E I ^ 0 (r 2 )7 10 (Q 1 ) >, (14) 

and 

D =< <S>f{r x ,r 2 ) \H-E\ <S>f{r x ,r 2 )>= s x -E, (15) 

s x is the expectation value of //, 

< | H | cD ( ,° >= £l . 


(16) 


andO^ is normalized to unity. The variation with respect to C\, i.e., 


51 

SC, 


= 0 , 


( 17 ) 


implies 
B + 2C X D = 0 . 


(18) 


This gives 


c B _ 4<O^ ) (r 1 ,F 2 )|//-£|^ 0 (r 2 )T 10 (Q 1 )> _ 2 <V x (r x )u{r x )> 
1 2D H-s j E-e l 


(19) 


S I 

The equation for the scattering function can be obtained from — — = 0 , which reduces to 

S u 


<Y l0 (Q l )</, 0 (5 2 )\H-E\W L >=0, (20) 

and this can be simplified to 

< ^ioC^iMjC^) \ H - E \ [(/> a (r 2 )u(r l )± (1 <-> 2)] > +C l < Y 10 (Q. l )^ 0 (r 2 ) \ H - E \ >= 0, 

(21) 

or 

< I H ~ E I [(^o(hMf)± (1 <-> 2)] + CjFjCrj) = 0, (22) 

Where E, (r, ) is given in Eq. (14). Substitution of Ci from Eq.(19) into Eq.(22) gives the 
equation 

2V (r)<V ( r)u( > 

< 7 10 (Q, M, (r 2 ) | H - E | 0 O (r 2 >#(/■ )±(lo2)]>+ ^ ^ = 0 (23) 

E-s j 


The resulting equation for « £ (r) can be written in the form 

( ^_L(i + l) + ^ ± ^_ F + ^ K(r) = 0> 

dr" r 


(24) 


where 


( 25 ) 


V l , L =2rY Vi(r)<Vl( '' l)U(r ' )> . 

°p L 77 
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In [6, 7], the optical potential of the type given in Eq. (24) was derived by using 
the Feshbach projection operator formalism [8] based on projection operators P and Q. 
The present formalism is independent of the projection operators P and Q. 

In the above fonnalism [6, 7], those terms coming from the correlation function 
are such that they take into account only the short-range correlations and there is no 
direct long-range potential proportional to 1/r 4 in the scattering equation satisfied by 

u{ r). 


III. Optical potential with polarization 

This long-range potential is due to the polarization of the target wave produced by 
the presence of the incident electron. This can be taken into account by the method of 
polarized orbitals. Temkin [12] has shown, using the adiabatic approximation in the first- 
order theory and using the dipole part of the resulting perturbed wave function, that in the 
presence of the incident electron r\ the effective target wave function can be written as 


® pol (r l ,r 2 ) = <j> 0 (r 2 ) 


ZpO \) Un^ P (r 2 ) cos($ n ) 
K r i 


where 
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+ r 2 ). 


(26) 


(27) 


andi9 12 is the angle between T\ and r 2 . We have replaced the step function s(r x , r 2 ) 
used by Temkin [12] by a smooth cutoff function X p( r \) which is of the form 


Xp( r \) = (\~e~P r ' )", (28) 

where n > 3. Now the polarization takes place whether the scattered electron is inside or 
outside the orbital electron. The polarization function given in Eq. (26) is valid 
throughout the range. This is unlike the step function s{r x ,r 2 ) used by Temkin [12] 
which ensures that the polarization takes place when the scattered electron r\ is outside 
the orbital electron r 2 . Furthermore, the function in Eq. (28) gives us another nonlinear 
parameter [], which is a function of k, the incident electron momentum, along with the 
exponent n. This term guarantees that X {S r \ V r \ 6 when r, — > 0 and it also 

contributes to the short-range correlations in addition to those obtained from the 


correlation function O l , and therefore, is useful to optimize the results. Now the electron- 
target wave function can be written as 




(29) 


We arrive at the same form of the scattering equation as Eq. (24), when we replace 
<f> 0 (r 2 ) by O ,ml (r ] , r 2 ) given in Eq. (26). We further restrict ourselves to L= 1 and we can 
write the final equation in the fonn 


+ k‘ + -4 + V, + v* ± (K+ FT )- v * Mn = 0 . 


P°1 \ T/ P°l ' 


dr 


.2 d 


op 


(30) 


We give below various quantities: 


D(r) = 


, 43 XpY) 2 

8 Z 6 - 2 ’ 


(31) 


The direct potentials are given by 


V d = 


Mzi) + 2 e - (z+ i), 

r r 


(32) 


and 


v d° l = Oi +x 3 ) + x 2 ^-, 


d 

dr 


(33) 


We give xi, X2 and X3 below 


Xj = 2 


(ZrY 


-a(r) 


(34) 


“(') = ! -f ^ + 7^) 3 + 13 W + 13.5(Z,) + 6.75) 


(35) 


& f \ 43 ,Xp(r) 
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*3 = + 43(5 2 (r) + T + k z ) 


Xp( r ) „ .Xp(r) 
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B l (r) = (nPe~ Pr (1 - e~ Pr ) n ~ l 


(37) 


2 Xp{r) 

— )/r 

r 


2 


B 2 (r) = n/3e Pr (. Bx(r ) + 5y(r)) / r 2 + 2( Xp [ ^ - B ^-) 

r r 

Bx{r ) = -(j8 + 2/ r)( 1 - 1 r ) ' 


43 48 

d(r) = + --e 2Z, '(0.5(Zr) 4 +4.5(Zr) 3 +18(Zr) 2 +45(Zr) + 81.25 + ^/xx) 

4Zr (Zr) 


dxx = 


106.75 96 48 

(Zr) + (Zr7 + (Z^ 


(38) 

(39) 

(40) 

(41) 

(42) 


It can be seen that (xi+ X3) has a term 


9 

2ZV ’ 


(43) 


where 9/(2Z 4 ) is the dipole polarizability of the target with nuclear charge Z. In addition 
to this term, there are other terms proportional to 1/r 2 , 1/r 3 , and 1/r 4 in the direct potential. 
They go to zero for r — >0. The exchange terms are given by 


V ex u(r) = 


8 Z 3 
3 


-Zr 


1 r 

[— [ dxe Zx x 2 u(x) + r 2 

V ^ 

r 0 


00 -Zx 

J dx — — u(x)]. 


(44) 


The exchange polarization terms are too detailed and are not given here. The 
optical potential is given by 




pol 


ft,Z 2 )\H-E\'¥ 0 > 


E-s „ 


(45) 


where 'Fq is the wave function V F 0 given in Eq. (29) without the correlation term O L . 

This Optical potential includes the effects of polarization. The right hand side of Eq. (45) 
has not been given explicitly because it contains a large number of terms. 


IV. Calculation and Results 


Eq. (30) is solved for the continuum function ui by the noniterative method of 
Omidvar [13] and the phase shift q is obtained from the value of the function at large 
distance: 


lim,^ u, (r) cc sin (hr -l — + rj). (46) 

In order to obtain phase shifts which can be compared to those obtained by the 
method of polarized orbitals, we exclude the correlation terms in Eq. (1). Also, the 
present cutoff function Zp given in Eq. (36), which allows optimization with respect to P 

to get the best results, is replaced by the cutoff function obtained by Shertzer and Temkin 
[14]: 

Zsr =l- e ( - 2Zr) [^ + ^p^ + 2(Zr) 2 +2(Zr) + l] (47) 

Now the calculation should be similar to that carried out by the method of 
polarized orbitals. The results obtained by the use of two different cutoffs are not very 
different. The phase shifts obtained, using this cutoff of Shertzer and Temkin [14], for l P 
and 3 P are given in Table I. and compared with those obtained by Sloan [15]. These 
results include contributions from the exchange polarization terms. The present results 
have variational bounds, i.e., they are always lower than the exact phase shifts. We see 
that phase shifts, obtained using the polarized orbital method, are always higher than the 
present ones. 

In Table II, we give the convergence of l P phase shifts for k=0 . 1 with increasing 
number of terms in the correlation function. We use the cutoff function given in Eq. (28). 
The results have been optimized with respect to the nonlinear parameters P, y and 5, with 
n= 3, the optimum value. We see that we do not need more than 35 terms to get results 
comparable to those obtained in [6] without the use of the polarization term in the target 
wave function and using the Feshbach formalism of projection operators, where it is not 
possible to modify the formulation in such a way as to produce a direct polarization 
potential proportional to 1/r 4 in the scattering equation. 

Gailitis [ 1 6] has shown that as the number of terms in the correlation function is 
increased the optical potential become more attractive. Consequently, phase shifts 
increase as the number of terms is increased. We see from the results given in Table II 
that this holds true. 

In Table III, we give l P phase shifts for values of the incident momentum from 
k= 0.1 to 0.8 for 35 terms and compare them with those obtained in [6] with 220 terms in 
the correlation function, but without the polarization term. We see that in most case the 
results are comparable in accuracy to those obtained previously [6] with longer 
expansions in the correlation function. This indicates that very long expansions do give 


fairly accurate results. We also compare the present results with the variational results of 
Ajmera and Chung [17], R -matrix results of Scholz et al. [4] and the finite element results 
of Botero and Shertzer [5]. Most of the results are comparable but the results obtained in 
[4, 5, 16] for k=0.4 to 0.8 are rather too low compared to the present results. The phase 
shifts obtained in the above-mentioned calculations do not have any variational bounds. 

It should be noted that the phase shifts increase up to k=0.3 and then start decreasing up 
to £=0.7 and increase again. 

In Table IV, we give the convergence of P phase shifts for £= 0.1 with the 
number of terms in the correlation function. Again, we use the cutoff function given in 
Eq. (28). The results have been optimized with respect to the nonlinear parameters [1, y 
and 5, with n= 3, the optimum value. We again see that we do not need more than 35 
terms to get results comparable in accuracy to those obtained in [6] without the use of the 
polarization tenn in the target wave function and using the Feshbach formalism of 
projection operators. 

In Table V, we give results P phase shifts for values of the incident momentum 
from £= 0. 1 to 0.8 for 35 terms and compare them with those obtained in [6] with 220 
terms in the correlation function, but without the polarization tern. We see that in most 
case the results are comparable in accuracy to those obtained previously [6] with longer 
expansions. This indicates that very long expansions do give fairly accurate results 
because of completeness in the whole space. We also compare the present results with the 
Kohn variational results of Armstead [18], R -matrix results of Scholz et al. [4] and the 
finite element results of Botero and Shertzer [5]. Most of the present results are 
comparable with the results obtained [4, 5, 14]. It should be noted that the phase shifts for 
3 P increase continuously throughout the range from £=0.1 to £=0.8. Although the 
differences are small, the present ’ P are the largest of the lower bound results and 
therefore are the best. 


V. Low energy scattering 

It is known [19] that at low energies L= 1 scattering, the long-range correlations 
contribute most to the phase shift: 

tan(/ 7 )/£ 2 = Ka/\5-Ak. (48) 

So that 

(tan(/ 7 (£j )/ £ 2 - tan(/ 7 (£ 2 ) / k; ) = -A(£, - £ 2 ) . (49) 

The first tenn in Eq. (48) is due to the long-range potential and the second term has 
contributions from the short-range correlations of which there are contributions from the 
cross of <f > po! and O i terms. Thus the values of A of the present calculation do not 
coincide with values obtained in our previous calculation [6] which included only the 
short-range correlations. Using phase shifts given in Table VI, we find 


(50) 


A t =- 1.002 and zl s = 2.942. 

O’Malley et al. [18] obtained ^(triplet) «- 1.3 and ^(singlet) « 1 .6 by fitting the phase 
shifts of the original polarized orbitals calculation (e-H) of Temkin and Lamkin [20]. But 
the latter calculation, although of historical importance, did not give precision results. 
Thus, I believe, the A values obtained in the present calculation are more accurate. 


VI. Conclusions 

In conclusion, we have applied the hybrid theory, in the presence of an optical potential, 
in which the long-range and short-range correlations, as in Eq. (29), can be taken into 
account at the same time. The close-coupling formalism with short-range correlations is 
like the present formalism, in the sense that the P states give the polarizability of the 
lower-lying 5-state of the target. The E-wave phase shifts are much more sensitive to 
polarization and short-range correlations than 5-wave phase shifts [7]. The present results 
are calculated variationally and therefore have lower bounds to the exact phase shifts. 

The present results are very close to those obtained in [6], using the Feshbach formalism 
in which it was not possible to define the projection operators P and Q to modify the 
target function to include the effect of the distortion produced by the incident electron. 
Moreover, shorter expansions of the order of 35 terms are needed to get results 
comparable in accuracy to those obtained in [6]. 
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Table I. Comparison of phase shifts r| without correlations with those obtained in [16] 



X P 

3 

P 

k 

n 

Ipo 

7 

Hpo 

0.1 

0.0057 

0.0067 

0.0094 

0.0109 

0.2 

0.0110 

0.0171 

0.0384 

0.0486 

0.3 

0.0006 

0.0210 

0.0867 

0.1151 

0.4 

-0.0090 

0.0163 

0.0148 

0.2005 

0.5 

-0.0295 

0.0064 

0.2100 

0.2867 

0.6 

-0.0495 

-0.0039 

0.2625 

0.3574 

0.7 

-0.0646 

-0.0100 

0.2999 

0.4063 

0.8 

-0.0721 

-0.0095 

0.3225 

0.4351 

[16] I. H. Sloan, 

Proc. Roy. Soc. 281 , 151 (1964). 


Table II. Convergence of l P phase shifts p for k= 0.1 with the number of terms. 


N 

B 

T 

6 

0 

0 

0.47 



0.0060892 

4 

0.27 

0.42 

0.99 

0.0062066 

10 

0.21 

0.30 

0.99 

0.0063271 

20 

0.34 

0.33 

0.84 

0.0063444 

35 

0.34 

0.30 

0.87 

0.0063508 


Table III. Singlet P phase shifts p for various k for A=35 


k 

H 

n a PQ 

n b AC 

Vrm 

Hfem 

0.1 

0.00635076 

0.0063083 

0.005782 

0.006 

0.006 

0.2 

0.01506556 

0.014988 

0.01445 

0.015 

0.0148 

0.3 

0.01670634 

0.016613 

0.01550 

0.016 

0.0160 

0.4 

0.01015347 

0.0099980 

0.00846 

0.009 

0.0090 

0.5 

-0.00061223 

-0.00084017 

-0.00287 

-0.002 

-0.0020 

0.6 

-0.01009367 

-0.010359 

-0.013029 

-0.012 

-0.0117 

0.7 

-0.01321557 

-0.013483 

-0.017225 

-0.016 

-0.0149 

0.8 

-0.00490138 

-0.0048524 

-0.009544 


-0.0068 

“Phase shifts obtained using t 

le Feshbach formalism [8]. 


b Variational results of Ajmera and Chung [17]. 

VCmatrix results of Scholz et al. [4]- 

d Finite element method results of Betero and Shertzer [5]. 


Table IV. Convergence of triplet P 

}hase shifts r| for k= 0.1 with the number of terms. 

N 

B 

T 

5 

Phase shift 

0 

0.45 

0.35 

0.80 

0.0098336 

4 

0.45 

0.35 

0.80 

0.0099701 

10 

0.45 

0.35 

0.85 

0.0103265 

20 

0.45 

0.37 

0.85 

0.0103714 

35 

0.38 

0.30 

0.80 

0.0103813 


Table V. Triplet P phase shifts r| for various k for N=2>5 


k 

H 

n a PQ 

n b Kh 

Wrm 

Wfem 

0.1 

0.01038234 6 

0.010382 

0.0101 

0.010 

0.0100 

0.2 

0.04536735 

0.045345 

0.0448 

0.045 

0.0452 

0.3 

0.1069312 

0.10679 

0.1059 

0.107 

0.1067 

0.4 

0.1888873 

0.18730 

0.1866 

0.187 

0.1873 

0.5 

0.2709762 

0.27058 

0.2700 

0.270 

0.2708 

0.6 

0.3416749 

0.34128 

0.3405 

0.341 

0.3417 

0.7 

0.3932100 

0.39257 

0.3918 

0.392 

0.3933 

0.8 

0.4277296 6 

0.42730 

0.425 


0.4283 

a Phase shifts obtained using t 

le Feshbach formalism [8]. 


b Kohn variational results obtained by Armstead [18]. 
e 7?-matrix results of Scholz et al. [4]- 
d Finite element method results of Botero and Shertzer [5]. 
e Phase shifts for k= 0. 1 and 0.8 are for 56 terms. 


Table VI. Low energy P- wave phase shifts for N= 35 and values of A from the effective 
range theory. 


State 

k\=0.04 

£ 2 =0.05 

A 

'P 

0.001303692 

0.001963464 

2.942 

6 P 

0.001564346 

0.002469346 

-1.002 



